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Chapter  1 

Introduction 


This  report  describes  the  2008  version  of  the  Atmospheric  Laser  and  Infrared 
Transmission  (ALIRT)  model  developed  by  LMI  for  the  Tactical  Air  Division  of 
the  Office  of  the  Secretary  of  Defense,  Program  Analysis  and  Evaluation.  ALIRT 
is  a  physical  optics  model  that  implements  a  method  for  calculating  the  peak1  and 
average  fluence  delivered  to  an  airborne  or  ground-based  target  from  a  laser 
weapon  on  an  airborne  platform,  when  the  beam  propagates  through  the  inhomo¬ 
geneous,  absorbing,  scattering  atmosphere  (refraction  and  extinction),  which  is  in 
turbulent  motion  (turbulence),  and  when  heating  of  the  air  by  the  beam  (blooming) 
may  have  significant  effects.  The  purpose  of  the  ALIRT  model  is  to  support  analy¬ 
sis  of  the  use  of  high-energy  lasers  and  infrared  equipment  in  air  combat  engage¬ 
ment  and  campaign  analyses. 

ALIRT  models  both  1.3152  micron  (p)  and  1.060p  wavelength  lasers.  There  is 
continuing  interest  in  the  potential  for  use  of  airborne,  tactical,  high  energy  laser 
weapons.  Early  efforts  in  this  area  focused  on  10. 6p  wavelength  carbon  dioxide 
(C02)  and  3.8  p  wavelength  deuterium  fluoride  (DF)  chemical  lasers.  Current  in¬ 
terest  is  focused  on  1.3152  p  chemical  oxygen  iodine  (COIL),  and  electrically- 
powered  solid  state  lasers.  The  most  mature  high  energy  solid  state  laser  technol¬ 
ogy  is  flashlamp-pumped  1.06  p  neodymium/glass  (Nd:glass),  however,  diode 
pumping,  and  alternative  lasing  molecules  and  substrates  are  under  active  consid¬ 
eration. 

ALIRT  user  inputs  include  laser  power,  telescope  diameter,  laser  platform  altitude 
and  speed,  target  bearing,  range  and  relative  speed,  target  incidence  angle,  and  a 
choice  of  atmospheres.  Outputs  include  target  height  from  a  spherical  earth;  dif¬ 
fraction-limited,  turbulence-scattered,  absorbed,  and  thermally  bloomed  intensity 
plots;  and  numerical  intensity  data.  For  modeling  simplicity,  we  maintain  separate 
ALIRT  models  for  1.3152  p  wavelength  iodine  and  1.060  p  wavelength  neodym¬ 
ium/glass  lasers. 

In  Chapter  2,  we  describe  ALIRT ’s  Microsoft  Excel  spreadsheet  and  Visual  Basic 
(VB)  structure.  Chapter  3  is  an  ALIRT  user  guide.  The  appendix  details  the  calcu¬ 
lations  of  peak  fluence  and  target  interaction. 


1  Here,  the  term  “peak”  refers  to  the  maximum  of  the  ensemble  average  fluence.  Peak  does 
not  refer  to  the  short-term  spot  peak  fluence.  Turbulence  causes  movement  of  the  beam,  which  can 
manifest  itself  either  as  smearing  (if  the  movement  is  small  relative  to  the  short-term  beam  diame¬ 
ter)  or  by  hopping  around  (if  the  movement  is  large  relative  to  the  short-term  beam  diameter). 
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Chapter  2 

ALIRT  Model  Structure  and  Flow 


ALIRT  is  implemented  in  an  Excel  workbook  containing  several  spreadsheets, 
and  the  model  includes  both  spreadsheet  and  VB  calculations.  This  chapter  de¬ 
scribes  the  spreadsheets,  VB  subroutines,  and  model  flow.  In  a  major  revision 
from  earlier  versions,  we  have  moved  many  large-array  calculations  from  the 
spreadsheets  into  VB,  significantly  increasing  execution  speed. 

ALIRT  Spreadsheet  Structure 


ALIRT  includes  a  primary  user  input/output  spreadsheet — “Main” — plus  spread¬ 
sheets  for  data,  calculations,  and  intermediate  results.  The  spreadsheets  also  in¬ 
clude  calculations  of  final  outcome  data  that  are  linked  to  the  “Main”  user 
input/output  spreadsheet.  Table  2-1  lists  and  briefly  describes  the  ALIRT  spread¬ 
sheets. 


Table  2-1.  ALIRT  Spreadsheets 


Name 

Description 

Main 

Primary  program  control  and  user  input/output  display 

test 

Intermediate  plot  data  output  and  plots 

wind 

Cross-beam  wind  calculation 

fulltraj 

Parallel  calculations  of  turbulence  effects  for  checks 

cn2 

Comparison  of  CnA2  models 

SR 

Unused — early  turbulence  calculations 

blooming 

Blooming 

xtn 

Calculation  of  extinction  for  0  to  15  km  altitude 

abs 

Calculation  of  molecular  and  aerosol  absorption  for  0  to  15  km  altitude 

attl  31 5_1 

1976  Standard  Atmosphere  molecular  and  aerosol  absorption  data 

att1315_1_m 

1976  Standard  Atmosphere  molecular  extinction  data 

attl  31 5 _ 2 

PcLnWin  Mid-latitude  Summer  molecular  and  aerosol  absorption  data 

att1315_2_m 

PcLnWin  Mid-latitude  Summer  molecular  extinction  data 

attl  3 1 5_3 

PcLnWin  Mid-latitude  Winter  molecular  and  aerosol  absorption  data 

att1315_3_m 

PcLnWin  Mid-latitude  Winter  molecular  extinction  data 

Please  note  that  the  generic  spreadsheet  name  “attl3 15  ”  is  used  for  both  1.3  p 
and  1.06  p  data. 
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ALIRT  Visual  Basic  Structure 

Primary  Subroutines 

Four  buttons  on  the  “Main”  spreadsheet  control  the  VB  code  in  ALIRT.  Three 
buttons  control  selection  of  the  modeling  atmosphere.  The  current  three  options 
are  1976  Standard  Atmosphere,  PcLnWin  Mid-latitude  Summer,  and  PcLnWin 
Mid-latitude  Winter.  Adding  atmospheres  would  be  straightforward.  The  fourth 
button,  “Geom,”  initiates  a  run  of  the  main  program,  by  calling  the  subroutine 
“runplots.”  Runplots  calls  the  subroutines  shown  in  Table  2-2. 

Table  2-2.  Primary  Subroutines  Called  by  Subroutine  Runplots 


Subroutine 

Description 

init 

Retrieves  input  parameters  for  laser  height,  target  range,  and  target  bearing. 

slickray 

Computes  refracted  ray  from  laser  to  target  by  integrating  the  beam  equations, 
starting  with  the  user-input  weapon  height  and  look  angle,  filling  the  arrays  x(i) 
and  h(i)  that  give  distance  along  the  refracted  ray  and  altitude  of  points  on  the 
ray. 

integral 

1  5 

jc^(h(zR)(l-z)Mz 

Approximates  the  turbulence  integral  z=0  numerically  using 

the  trapezoidal  rule.  The  subroutine  assumes  x(i)  and  h(i)  are  in  place,  from 
slickray.  Results  are  placed  in  global  variable  “integ  “for  use  by  other  routines. 

plotfreespace 

Plots  the  intensity  in  kilowatts/cmA2  of  the  TEMOO  Gaussian  beam  propagat¬ 
ing  in  free  space,  and  the  same  quantity  for  that  beam  propagating  through 
turbulence  and  impacted  by  atmospheric  extinction  and  thermal  blooming,  at 
normal  incidence  on  the  target.  The  subroutine  uses  the  value  of  integ,  pro¬ 
duced  by  subroutine  “integral.”  Extinction  and  blooming  data  are  taken  from 
sheets  “xtn”  and  “bloom,”  respectively.  Peak  axis  intensity  values  for  a  beam 
affected  by  diffraction,  by  turbulence,  by  extinction,  and  by  blooming  striking  a 
target  at  normal  incidence  are  output  to  “Main.”  Plot  contour  values  for  these 
intensities  and  for  blooming  are  output  to  “test”  and  plotted,  and  the  plot  is 
copied  to  “Main.” 

plotcontours 

Plots  three  contours  of  constant  intensity,  for  a  TEMOO  beam  affected  by  tur¬ 
bulence,  extinction,  and  blooming  incident  at  a  user-specified  incident  angle 
on  a  plane  surface.  The  contours  are  output  to  “test,”  and  the  plot  is  copied  to 
“Main.”  The  peak  intensity  value  is  output  to  “Main.”  The  peak  intensity  for  the 
bloomed  beam  at  normal  incidence  is  output  to  “Main”  from  “plotcontours.” 

The  subroutine  assumes  public  variables  wturb  and  paxisturb  have  been 
evaluated  by  subroutine  plotfreespace. 
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ALIRT  Model  Structure  and  Flow 


Primary  Functions 

In  addition  to  the  subroutine  structure,  ALIRT  uses  several  VB  functions  in  the 
calculations.  Table  2-3  lists  and  describes  them. 


Table  2-3.  Primary  Functions  Used  by  ALIRT 


Function 

Description 

stda_ndx 

Refractive  index  of  1976  Standard  Atmosphere 

nmwl_ndx 

Refractive  index  of  Mid-latitude  Winter  Atmosphere 

xmls_ndx 

Refractive  index  of  Mid-latitude  Summer  Atmosphere 

dstda_ndx 

Derivative  of  Standard  Atmosphere  refractive  index 

dnmwl_ndx 

Derivative  of  Mid-latitude  Winter  refractive  index 

dxmls_ndx 

Derivative  of  Mid-latitude  Summer  refractive  index 

Huff_vall57 

C2 

Huffnagel-Valley  model  of  n 

clearl 

C2 

Clear  1  night  model  of  n 

ALIRT  Calculation  Flow 


On  “Main”  spreadsheet,  the  user  chooses  an  atmosphere  option  and  enters  laser 
power  and  mirror  diameter.  The  user  specifies  the  engagement  geometry,  giving 
the  weapon  height,  the  range  to  the  target,  the  target’s  bearing  and  elevation,  and 
the  angle  of  inclination  of  the  laser  beam  on  the  target  surface.  It  is  important  to 
note  that  the  target  elevation  angle  and  range  inputs  are  the  elevation  and  range  of 
the  refracted  ray  seen  at  the  laser;  they  are  not  the  elevation  angle  or  the  range  of 
the  straight  line  to  the  target. 

The  user  also  inputs  the  platform’s  airspeed  and  the  target’s  velocity  relative  to 
the  platfonn.  Relative  velocity  inputs  are  coordinates  of  target  velocity  in  a  rec¬ 
tangular  Cartesian  reference  frame  moving  with  the  weapon  aircraft,  which  is  as¬ 
sumed  to  be  in  straight,  level,  unaccelerated  flight.  The  “axial,”  “crossing,”  and 
“vertical”  axes  coincide,  respectively,  with  the  weapon  aircraft’s  roll,  pitch,  and 
yaw  axes.  All  axes  are  positive  away  from  the  aircraft. 

When  all  inputs  are  made,  the  user  clicks  “Geom,”  which  invokes  a  sequence  of 
VB  subroutines  and  spreadsheet  calculations.  The  first  of  these  is  a  calculation  of 
the  refracted  ray’s  path  through  the  inhomogeneous  atmosphere.  This  is  done  by 
numerical  integration  of  the  ray  equation,  following  a  method  of  Roland  Berton. 1 
We  neglect  earth’s  oblateness  and  assume  that  refractive  index  is  a  function  of 
height  only.  We  do  not  consider  earth-grazing  effects. 

1  Roland  P.H.  Berton,  “Variational  Calculation  of  Three-Dimensional  Atmospheric  Refrac¬ 
tion:  Part  I,  Description  and  Validation  of  the  Method,”  Journal  of  Optics  A:  Pure  and  Applied 
Optics,  August  2006,  pp.  817-830. 
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We  use  the  refracted  ray  rather  than  a  straight  line  in  space  for  two  reasons.  First, 

C2 

refractive  effects  can  cause  noticeable  variations  in  the  turbulence  parameter  n 
even  though  they  cause  only  small  departures  from  the  straight  line,  because 
r2 

n  varies  strongly  with  height  at  certain  altitudes.  Second,  we  want  to  be  consis¬ 
tent  with  Ontar  Corporation’s  PcLnWin  calculations,  which  use  refracted  paths. 

With  the  refracted  ray  path  computed,  ALIRT  calculates  the  diffraction-limited 
intensity  pattern  at  the  target.  We  use  the  standard  paraxial  approximation  for 
free-space  propagation  of  a  TEMOO  Gaussian  beam,  which  results  in  a  Gaussian 
beam  on  the  target.2  We  realize  that  the  assumption  of  a  TEMOO  beam  is  a  heavy 
one,  in  view  of  the  many  features  of  real  lasers  that  militate  against  such  a  simple 
beam.  Nevertheless,  we  make  that  assumption,  deferring  introduction  of  more 
complex  beam  models  until  we  have  data  on  the  beam  structures  of  specific  laser 
weapons. 

Next,  ALIRT  evaluates  the  impact  of  atmospheric  turbulence  on  the  beam.  The 
calculations  of  turbulence  effects  use  methods  given  in  The  Infrared  Handbook ,3 
The  calculations  involve  an  integral  along  the  center  axis  of  the  refracted  beam, 
which  ALIRT  approximates  by  the  trapezoidal  rule.  The  end  result  is  a  prediction 
for  ensemble  averages  of  the  turbulence-impacted  beam  intensity  on  the  target  as 
Gaussian,  with  a  larger  diameter  than  that  of  the  free-space  TEMOO  beam. 

ALIRT  then  evaluates  the  effects  of  atmospheric  extinction.  We  do  this  by  devel¬ 
oping  the  atmospheric  extinction  coefficient  as  a  function  of  height  above  a 
spherical  earth  and  then  integrating  along  the  refracted  beam  axis.  In  essence,  the 
atmospheric  extinction  function  is  an  interpolation  in  tabular  data  for  the  loga¬ 
rithm  of  the  extinction  coefficient,  which  we  obtain  from  Ontar’ s  PcLnWin  appli¬ 
cation  package.4  The  tabular  data  from  PcLnWin  are  contained  in  spreadsheets 
attl 3 1 51,  attl315_lm,  att  131 5_2,  attl315_2m,  attl 3 1 5  3,  and  attl315_3m.  The 
index  1  indicates  PcLnWin  data  for  the  1976  Standard  Atmosphere,  index  2  de¬ 
notes  data  for  Mid-latitude  Winter,  and  index  3  indicates  data  for  Mid-latitude 
Summer.  The  index  m  denotes  molecular  absorption  data,  used  in  blooming  cal¬ 
culations.  We  assume  that  extinction  attenuates  the  beam  without  distorting  the 
shape  of  the  ensemble-averaged  intensity  variation. 

Finally,  ALIRT  models  the  impact  of  blooming.  The  calculations  follow  the 
method  of  Weichel,  as  given  in  the  American  Physical  Society’s  report  on 
boost-phase  intercept  systems  for  national  missile  defense.5  Energy  input  for 


2  The  calculations  are  described  in  LMI,  2007  Atmospheric  Laser  and  Infrared  Tracking 
Model:  Some  Considerations  for  Interactions  of High-Power  Laser  Beams  with  Targets ,  Report 
PA603T3,  David  A.  Lee,  September  2007;  relevant  portions  are  in  the  appendix  to  this  document. 

3  See  Note  2. 

4  The  modeling  of  the  atmospheric  extinction  coefficient  is  detailed  in  LMI  Report  PA603T3, 
relevant  portions  of  which  are  in  the  appendix  to  this  document. 

5  American  Physical  Society,  Report  of  the  American  Physical  Society  on  Boost-Phase  Inter¬ 
cept  Systems  for  National  Missile  Defense,  July  2003. 
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ALIRT  Model  Structure  and  Flow 


blooming  is  taken  to  be  molecular  absorption,  obtained  from  PcLnWin  and  pro¬ 
vided  in  spreadsheets  attl315_xm,  where  x  is  1,  2,  or  3,  as  explained  above.6 

Our  models  have  certain  limitations.  We  neglect  earth’s  oblateness  and  do  not 
consider  earth-grazing  effects  on  propagation.  Our  model  of  atmospheric  extinc¬ 
tion  extends  only  to  1 5  km  height.  ALIRT  was  originally  developed  to  study  laser 
weapons  for  tactical  air-to-air  engagements,  for  which  the  15  km  limit  is  ade¬ 
quate. 

A  box  on  the  “Main”  spreadsheet  shows  the  beam’s  minimum  altitude;  the  user  is 
thus  warned  if  grazing  might  be  a  problem.  The  beam’s  maximum  altitude  also  is 
displayed.  ALIRT  may  be  used  for  cases  involving  greater  heights,  but  the  user 
should  note  that  the  results  will  be  those  for  which  atmospheric  extinction  disap¬ 
pears  at  altitudes  above  15  km. 


6  Further  details  on  the  blooming  calculations  are  given  in  the  appendix. 
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Chapter  3 

ALIRT  User  Guide 


All  ALIRT  operations  use  the  “Main”  spreadsheet.  Figure  3-1  shows  the  spread¬ 
sheet’s  control  buttons,  input  numbers,  and  output  numbers,  and  Figures  3-2  and 
3-3  show  the  output  plots. 

Figure  3-1.  “Main"  Spreadsheet:  Control  Buttons,  Input  Numbers, 
and  Output  Numbers 


PcLnWin  Mid-Latitude  Summer 

In  use: 
1976  U.  S. 

Standard 

Inputs 

Laser  Properties 

Engagement  Geometry  Geom  | 

Speeds,  kt 

Power,  MW 

0.5 

Weapon  height,  km 

5 

Platform  400 

Mirror  Diameter,  m 

0.5 

Range,  km 

35 

Target  relative  speeds 

Bearing,  degrees 

90 

Rel.  axial  0 

Wavelength  p 

1.3152 

Elevation,  degrees 

0 

Rel  crossing  0 

(Currently  Iodine  only) 

Incidence  angle,  degrees 

PLEASE  CLICK  "Geom"  BUTTON  AFTER 
CHANGING  ANY  INPUT 

45 

Rel.  vert  0 

Elevation  angle  for  ground  target  at  entered  range,  unrefracted  beam 

-8.369  0 

Earth  tangent  elevation  angle  and  range,  unrefracted  beam 

-2.27  0 

Outputs 

Peak  intensity  on  target,  kw/(cmA2) 

Intensities  for  normal  incidence, 

for  angled  incidence,  with  turbulence, 

at  center  of  spot,  kw/(cmA2): 

Target  height,  km  Min.  beam 

extinction,  and  blooming. 

Diffraction  limited 

4.63E+00 

5.09  5.00 

With  atmospheric  turbulence 

2.90E+00 

Cautions: 

1.35E+00 

With  atmospheric 

1.  The  current  model  only  includes 

turbulence  and  extinction 

2.62E+00 

therefore,  results  may  be  optimi: 

With  turbulence,  extinction 

2.  Turbulence  below  1  Km  is  based 

and  blooming 

1.91E+00 

local  effest  of  terrain,  foliage,  et( 

3.  Results  unreliable  for  Earth-grazi 

4.  All  propagation  is  for  cloud-free, 
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Figure  3-2.  “Main”  Spreadsheet:  Normal  Incidence  Intensity  Profile  Plots 


Figure  3-3.  “Main"  Spreadsheet:  Selected  Incidence  Angle  Spot  Size  Plots 

Constant-intensity  contours 
for  Angled  Incidence 
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ALIRT  User  Guide 


Using  ALIRT  requires  the  following  actions: 

1 .  Set  atmosphere  by  clicking  the  button  beside  the  desired  atmosphere. 
ALIRT  posts  the  name  of  the  chosen  atmosphere. 

2.  Enter  laser  power  in  megawatts  and  mirror  diameter  in  meters. 

3.  Enter  weapon  height  and  range  in  kilometers. 

4.  Enter  target  bearing  in  degrees.  Bearing  is  the  angle  between  the  weapon 
aircraft’s  flight  direction  and  the  line  of  sight  to  the  target,  which  is  as¬ 
sumed  to  be  to  port. 

5.  Enter  target  elevation  in  degrees.  Elevation  is  the  angle  between  the  re¬ 
fracted  laser  ray  to  the  target  as  it  leaves  the  weapon  aircraft  and  the  plane 
as  determined  by  the  aircraft’s  roll  and  pitch  axes. 

6.  Enter  the  angle  of  incidence  of  the  laser  ray  on  the  target.  Perpendicular 
incidence  is  0  degree.  ALIRT’ s  calculations  do  not  treat  grazing  effects,  so 
the  incidence  angle  should  not  be  greater  than  60  degrees  in  magnitude. 

7.  Enter  the  weapon  aircraft’s  true  airspeed  in  knots. 

8.  Enter  the  components  of  the  target  aircraft’s  velocity  relative  to  the 
weapon  aircraft,  in  knots:  “axial”  component  is  in  the  direction  of  the 
weapon  aircraft’s  roll  axis,  “crossing”  component  is  in  the  direction  of  the 
weapon  aircraft’s  pitch  axis,  and  “vertical”  component  is  in  the  direction 
of  the  weapon  aircraft’s  yaw  axis. 

When  these  entries  are  complete,  click  “Geom.”  ALIRT  then  posts  these  outputs: 

1 .  Plots  of  intensity  vs.  radial  distance  from  the  center  of  the  beam  for  the 
diffraction-limited  beam  and  for  the  beam  affected  by  turbulence,  extinc¬ 
tion,  and  blooming. 

2.  Intensities  at  the  center  of  the  beam  at  normal  incidence  on  the  target,  in 
kw/cm2,  for  the  diffraction-limited  beam;  for  the  beam  affected  by  diffrac¬ 
tion  and  turbulence;  for  the  beam  affected  by  diffraction,  turbulence,  and 
extinction;  and  for  the  beam  affected  by  diffraction,  turbulence,  extinction, 
and  blooming. 

3.  Target  height  and  minimum  height  of  the  refracted  beam.  The  minimum 
height  can  be  less  than  either  weapon  height  or  target  height  for  ranges 
greater  than  the  earth-grazing  distance. 

4.  The  elevation  angle  of  the  straight  line  from  the  weapon  to  a  ground  tar¬ 
get.  The  elevation  angle  of  the  refracted  beam  striking  a  ground  target  will 
be  slightly  greater  (algebraically)  than  this  value.  For  the  data  shown  in 
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Figure  3-1,  the  earth-strike  elevation  of  the  refracted  ray  is  -5.855  de¬ 
grees,  rather  than  the  -5.873-degree  elevation  of  the  straight  line. 

5.  The  elevation  angle  and  length  of  the  straight  line  from  the  weapon  that  is 
tangent  to  the  earth.  This  information  is  intended  to  help  the  user  avoid 
grazing  beams,  which  ALIRT’s  calculations  do  not  treat. 
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Appendix 

Calculations  of  Peak  Fluence  and  Target 
Interaction 


This  appendix  consolidates  and  updates  material  from  LMI  reports  [1]  and  [2], 
Certain  sections  describe  significant  improvements  in  the  model. 

The  appendix  begins  with  a  brief  overview.  We  then  derive  the  field  equations  we 
use,  starting  with  Maxwell’s  equations  and  identifying  our  assumptions.  Next,  we 
develop  an  integral  equation  for  electromagnetic  propagation  in  a  random  medium 
occupying  a  half-space,  with  field  values  prescribed  on  the  bounding  plane.  The 
integral  equation  allows  us  to  describe  diffraction  effects  and  to  address  the 
effects  of  imperfect  initial  beam  quality.  It  also  serves  as  our  starting  point  for 
turbulence,  extinction,  blooming,  and  other  calculations.  A  mathematical  note  at 
end  of  the  appendix  develops  expressions  for  the  electric  field  and  for  intensity 
(power  density). 

Modeling  Overview 


ALIRT  treats  the  effects  of  beam  quality,  diffraction,  atmospheric  turbulence, 
atmospheric  absorption  and  scattering,  and  blooming  on  the  energy  delivered  at 
not-necessarily-normal  incidence  on  a  not-necessarily-plane  target. 

We  do  not  consider  improvements  by  adaptive  optics  for  two  reasons.  First,  we 
expect  the  tactical  engagements  of  interest  to  take  place  at  ranges  of  no  more  than 
a  few  tens  of  kilometers,  reducing  the  impacts  of  atmospheric  turbulence 
significantly  in  comparison  with  those  of  strategic  laser  weapon  engagements 
expected  to  occur  at  ranges  of  some  hundreds  of  kilometers.  Second,  we  want  to 
see  what  might  be  done  with  tactical  laser  weapons  without  the  complexity, 
weight,  initial  cost,  and  maintenance  expense  associated  with  adaptive  optics. 

The  extensive  literature  on  Gaussian  beams  treated  with  the  classical  Rytov  model 
of  turbulence-induced  fluctuations  in  refractive  index  allows  us  to  describe  the 
combined  effects  of  diffraction  and  turbulence. 

We  developed  models  of  extinction  coefficient  and  molecular  absorption  from 
PcLnWin  outputs  in  Ontar  [3].  We  used  these  models  to  calculate  effects  of 
atmospheric  extinction  and  blooming. 

Much  of  the  material  in  this  appendix  is  available  in  standard  references, 
particularly  Andrews  and  Phillips  [4]  and  the  references  cited  in  that  work.  For 
our  discussions  of  beam  quality  and  of  non-nonnal  incidence,  we  need  certain 
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details  of  Gaussian  beam  propagation.  For  completeness  and  clarity,  we  derive  the 
results  that  we  need.  This  makes  our  notation  conventions  clear  (no  wondering  if 
wave  number  is  the  reciprocal  of  wavelength  or  2n  times  that  value)  and  explicitly 
identifies  the  important  assumptions. 

The  Field  Equations 


In  a  medium  with  constant  permeability  p0  and  variable,  scalar  pennittivity  s,  and 
no  true  currents  or  true  charges,  Maxwell’s  equations  are 


V  •  (sE)  =  0 

V  -H  =  0 

VxE  =  -p0  — 
dt 

VxH  =  « 
dt 


[Eq.  A-l] 


As  is  usual  in  beam  propagation  analyses,  we  assume  that  spatial  gradients  and 
time  derivatives  of  s  are  negligibly  small  compared  with  those  of  the  field 
components.  We  also  set 

U(x1,x2,z,t)  =  U(x1,x2,z)e_Kot  [Eq.  A-2] 


where  U  is  either  of  E  or  H  . 

With  those  assumptions,  Equations  A-l  imply  that  both  electric  and  magnetic 
fields  are  solenoidal  and  that 


VxE  =  icop0H 


and 


V  x  H  =  -icosE  . 


[Eq.  A-3] 


[Eq.  A-4] 


Taking  the  curl  of  Equation  A-3,  using  the  fourth  of  Equations  A-l  to  express  the 
curl  of  H  in  terms  of  E ,  and  recalling  that  the  curl  of  the  curl  of  a  solenoidal 
vector  field  is  the  negative  of  the  Laplacian  of  the  field,  we  find  that  the  electric 
field  satisfies 


V2E  + co2p0sE  =  0  [Eq.  A-5] 

at  interior  points  of  the  medium. 


A-2 


Calculations  of  Peak  Fluence  and  Target  Interaction 


Introducing  a  spatially  varying  electric  susceptibility  y,  we  find 

V2E  + co2p0s0(l  + y)E  =  0  [Eq.  A-6] 


or 


V2E  +  k2(l  +  y)E  =  0 


[Eq.  A-7] 


where 


22  UJ 

k  =co  p0s0  =  — 
C0 


[Eq.  A-8] 


and  Co  is  the  speed  of  light  in  free  space. 

An  Integral  Equation  for  a  Field  Coordinate 


We  are  interested  in  electric  field  coordinates  U(xi,X2,z)  that  satisfy  Equation  A-7 
in  the  half-space  z  >  0,  with  the  boundary  condition 

U(x1,x2/0)  =  U0(x1,x2),  [Eq.  A-9] 

where  Uo(xi,  X2)  is  a  prescribed  boundary  value  function,  and  a  radiation 
condition  specifying  that  energy  propagated  by  the  electromagnetic  field  moves 
toward  increasing  values  of  z  for  sufficiently  large  z.  Such  a  field  coordinate  will 
satisfy  the  equation 


V2U  +  k2U  =  -k2yU  [Eq.A-10] 

and  the  boundary  condition  Equation  A-9. 

The  boundary  value  problem  posed  by  Equation  A-9  and  the  equation 

V2U  +  k2U  =  f(x!,x2,z)  [Eq.  A- 11] 

may  be  solved  with  a  Green’s  function  G(x,y)that  has  the  following  conditions: 

1 .  V2G  +  k2G  =  0  in  R+,  where  R+  denotes  the  half-space  y3  >  0  with  the 
point  y  =  x  excluded;  here  V“  denotes  the  Laplacian  operator  applied  to 
the  components  of  y  . 

e^-y\  _  ,  . 

2.  G(x,y)  =  — - ho(|x-y|  )  as  |  x  -y  |  T  0,  V  x  with  x3  >  0  . 

I x — y  I 
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3.  G(x,y)  s  0,  when  y3  =  0 . 

When  such  a  function  can  be  found,  the  solution  of  the  boundary  value  problem 
posed  by  Equation  A-9  and  Equation  A-l  1  may  be  written  as 

4rcU(x)  =  JjjG(x,y)f(y)dvy  -  JJ  ^^U0(y1,y2)dy1dy2  .  [Eq.A-12] 

R+  y3=0  vyi 

It  is  not  difficult  to  show  that  the  function 

Pikb-yl  Pikh-y-l 

G(x,y)  =  — - -  -  — - -  [Eq.  A- 1 3] 

|  x-y  |  |  x-y_  | 

satisfies  conditions  1,  2,  and  3,  where  y_  =  (yi/y2,-y3)  •  Consequently,  a  field 

component  U  that  satisfies  Equation  A-9  and  Equation  A- 10  satisfies  the  integral 
equation 


47rU(x)  =  -k2jjjG(x,y)x(y)U(y)dvy  - 

R+ 

JJ  ^^U0(yI/y2)dy1dy2 

y3=0  3 

Diffraction  and  Beam  Quality 


[Eq.  A- 14] 


We  will  study  diffraction  and  beam  quality  with  Equation  A- 14,  specialized  to  the 
case  x  =  0.  In  this  case,  Equation  A- 14  gives 


-I  oo  oo 

uoo=~n 


<3G(x,y) 


47r  _CxL  fy3 


u0(yi,y2)dy1dy2 . 


[Eq.  A- 15] 


Differentiating  the  function  G(x,y)  of  Equation  A- 13  with  respect  to  y3  gives 
<9G(x,y) 


dy3 


rik 

_o 

=  -2 

9 

y  3=0 

yro 

X3  cikrn 


[Eq.  A- 16] 


where 


r0  ^V(xi-yi)2  +  (x2-y2)2  +  x3  •  [Eq-  A-i7] 

We  are  interested  in  the  beam  far  from  the  initial  plane  x3  =  0.  We  assume  that  the 
beam’s  intensity  is  negligible  except  inside  a  right  circular  cone  of  angle  df 
whose  axis  is  the  x3  axis  and  where  VE  does  not  exceed  a  few  tens  of  microradians. 
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(The  following  analytic  results  can  be  checked  for  consistency  with  this 
assumption,  and  observations  support  its  relevance  to  actual  beams.) 

With  this  assumption,  for  all  values  of  xi,  yi  of  interest, 

(X‘~yi')  +(X2~y2)  «i.  [Eq.  A- 18] 

*3 

Using  this  assumption,  writing  r0  as 


r0  =x3Jl  + 


(xi  -yi)“  +(x2 -y2y 


[Eq.  A- 19] 


(x,  —  y,)2  +  (x?  -yd2 

expanding  the  radical  in  powers  of - - - - — - — — — ,  retaining  the  leading 

X3 

term  in  the  denominators  of  Equation  A- 16,  and  retaining  the  two  leading  terms  in 
the  exponential  function  of  that  equation, 1  we  find 


<9G(x,y) 


cry  3 


2ik 


-exp 


y3=o 


ik  ik 


+^— [(xi-yi)2+(x2-y2)2 


x3  2x3 


[Eq.  A-20] 


so  that 


U(x)  =  - 


ik 


ikx, 


2nx-, 


f  I  exp 


ik 

2x, 


(xi  -  yi)2  +(x2  -y2)2 


uo(y^y2)dyidy2  •  [Eq- A-21] 


Equation  A-21  gives  values  of  the  field  coordinate  U(x)  far  from  the  initial  plane 
and  allows  us  to  consider  the  effects  of  beam  quality  and  diffraction. 


Effects  of  Beam  Quality 


Under  certain  conditions  on  U0(xi,  x2),  which  we  state  below,  the  propagating 
field  will  have  relatively  simple  characteristics,  such  as  those  of  the  TEMoo  mode 
But  some  features  of  high-power  lasers  militate  against  this.  In  a  transverse-flow 
gas  laser,  upstream  depletion  of  the  inversion  may  destroy  axial  symmetry  of  the 
output.  Support  structures  for  optics  may  block  parts  of  the  output  beam. 
Vibrating  structures  (jitter)  may  spoil  intended  symmetries  or  distort  the  beam. 


1  The  additional  term  is  kept  in  the  exponent  because  kxpan'F  may  be  large  compared  with  1, 
and  so  affect  phase  in  the  beam,  while  kx3(tanvF)4  will  be  small  compared  with  1.  For  example,  if 
the  beam  is  negligible  for  distances  greater  than  1  meter  from  its  axis  at  a  range  of  20  km,  so  that 
*P  =  50  pradian,  when  X  =  1.3  p,  kx3  tan*F  =  O(106),  while  kx3(tan'F)4  =  o(10'7). 
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To  develop  some  understanding  of  departures  from  TEM0o,  we  suppose 
thatU(x)  is  the  coordinate  of  a  rectilinearly  polarized  electric  field,  and  that 


U0(x1/x2)  =  — 


i^n 


Pmn^O 


2m+n-27rm!n! 


H, 


V2x,  V2x2 

1(^77±)Hn(^T7^)exP 


wn 


wn 


r  1 

ik  'l 

2 

— 

— — 

r 

IWo 

2F0; 

[Eq.  A-22] 


where  the  sums  range  over  a  set  of  values  of  m  and  n,  where  r2  =  x2  +  x2 ,  where 

R0  =  I—  =  377  f2  is  the  “resistance  of  free  space,”  and  where  c))mn  is  a  phase 
V  £o 

factor  chosen,  with  Pmn,  to  make  the  sum  represent  the  initial  field.  The  Hj  are 
Hermite  polynomials  of  order  j. 


This  is  not  the  most  general  expression  for  the  field  boundary  values — the  widths 
in  the  xi  and  X2  directions  need  not  be  equal — but,  by  virtue  of  the  completeness 
property  of  Hermite  polynomials,  it  is  capable  of  representing  a  reasonably 
diverse  class  of  initial  fields.  The  components  of  the  sum  may  be  called 
“components”  of  the  input  field. 

Equation  A-22  represents  an  initial  field  that  has  been  modified  by  reflection  in  a 
thin  parabolic  mirror  with  focal  length  Fo.  By  virtue  of  the  orthogonality  property 
of  the  Hermite  polynomials,  the  total  power  P  in  the  input  is  the  sum  of  the 
powers  Pmn  of  the  individual  components. 

For  the  rectilinearly  polarized  electric  field  that  we  are  considering,  the  Poynting 
vector  will  point  along  the  X3-axis,  and  the  magnetic  field,  like  the  electric  field, 
will  be  perpendicular  to  that  axis.  It  is  customary  to  speak  of  the  component  with 
indices  m,  n  as  a  “Gauss-Hennite  TEMm,n  mode.” 

In  the  mathematical  note  at  the  end  of  this  appendix,  we  develop  expressions  for 
the  electric  field  and  for  intensity  (power  density)  at  distances  X3  from  the  initial 
plane,  for  the  Gauss-Hennite  TEMm  n  mode.  In  the  focal  plane  X3  =  F0,  the  field  is 


U(x1,x2,F0)  =  -|^^e“ 
w 


brnn^O 


2m+n~2  jim!n! 


V2 x .  V2 x , 

Hm(A-^)Hn(A-^)exp 


W 


W 


1 


ik 


W2  2F, 


oy 


[Eq.  A-23] 


where 


2F 

W  =  — .  [Eq.  A-24] 

kW0 

The  results  show  that  even  relatively  modest  diversions  of  power  from  the  TEM0,o 
mode  can  have  significant  impacts  on  the  peak  intensity  at  the  target  of  a  focused 
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beam.  Figure  A-l  and  Figure  A-2  give  an  arbitrary,  illustrative  example.  Figure 
A-l  shows  the  two-dimensional  intensity  profile  when  all  the  beam’s  power  is  in 
the  TEM0,o  mode,  and  Figure  A-2  gives  the  intensity  profile  when  10  percent  of 
the  power  is  in  the  T EM 2,2  mode.  The  peak  intensity  is  only  60  percent  as  great  in 
the  second  case  as  in  the  first. 

Figure  A-1.  Intensity  Profile  of  TEM0,o  Beam 


Figure  A-2.  Intensity  Profile  of  Beam  with  90  Percent 
of  Power  in  TEM0:o,10%  in  TEM2,2 


The  fact  that  modest  departures  from  a  pure  GaussianTEMo.o  mode  can  have 
significant  effects  should  not  be  surprising,  since  the  TEMo,o  mode  is  a  highly 
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idealized  case.  Pampaloni  and  Enderlein  [5]  have  shown  that  it  is  the  beam  with 
simultaneously  minimal  divergence  and  minimum  transverse  extension. 

Initial  conditions  other  than  Gaussian  beams  have  been  considered  in  the  laser 
propagation  literature.  For  example,  the  American  Physical  Society  (APS)  study 
of  boost-phase  intercept  systems  [6]  implicitly  assumes  that  Uo(xi,  X2)  is  a 
truncated  spherical  wave.  This  assumption  also  underlies  some  of  the  work  in 
LMI  Report  PA304T3. 

We  will  continue  our  work  assuming,  as  do  many  studies,  that  all  the  beam’s 
energy  is  in  the  TEM0jo  mode.  But  we  note  that  achieving  this  in  an  actual  device 
is  a  challenge. 

Effect  of  Diffraction  on  a  TEM0,o  Beam 

Figure  A-l  illustrates  the  effect  of  diffraction  on  a  TEMo.o  Gaussian  beam. 
Specifically,  as  shown  in  the  mathematical  note,  at  the  focal  plane  X3  =  F,  this 
beam’s  intensity  is  given  by 


I  = 


7tW 


where  Poo  is  the  total  power  in  the  beam,  and 


w  = 


V2F0  _  1  X 

kw7“ 


Fn . 


[Eq.  A-25] 


[Eq.  A-26] 


Relating  the  reference  length  Wo  of  the  Gaussian  beam  to  the  diameter  D  of  the 
mirror  that  launched  the  truncated  spherical  beam  considered  by  the  APS  so  as  to 
make  the  on-axis  intensity  the  same  in  both  cases  implies 


W 


0  ~ 


_D_ 

V8 


[Eq.  A-27] 


and  leads  to  the  comparison  between  intensity  patterns  of  the  Gaussian  and 
truncated  spherical  beams  shown  in  Figure  A-3.  The  relatively  close  agreement 
suggests  that  considering  Gaussian  TEMoo  beams  for  damage  effects  of  high- 
power  laser  beams  is  not  likely  to  lead  to  different  conclusions  from  considering 
truncated  spherical  beams,  as  was  done  in  LMI  Report  PA304T3. 
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Figure  A-3.  Intensity  Comparison 


- Gauss-  Tr.  Sph. 


Turbulence 


To  account  for  the  effects  of  atmospheric  turbulence,  we  return  to  Equation  A- 14 
written  in  the  fonn 


U(x)  =  -f-fflG(x,y)z(y)U(y)dvy  +Uf(x)  [Eq.  A-28] 

4;r  R+ 

where  Uf  is  the  free-space  solution  given  by  the  U  of  Equation  A-23. 

It  will  prove  convenient  to  write  Equation  A-28  in  the  fonn 

U(x)  =  H[U ]  +  Uf  (x)  [Eq.  A-29] 

where  the  operator  H  is  defined  by 

H[U(x)]  =  “^~JjjG(x,y)x(y)U(y)dvy  .  [Eq.  A-30] 

^71  T?  + 


Equation  A-28  is  a  second-kind  Fredholm  equation,  a  member  of  a  widely  studied 
class.  In  particular,  arguments  such  as  those  given  by  Oden  [7]  show  that, 
provided  the  largest  value  of  |y|  in  the  beam  region  is  sufficiently  small,  the 
sequence  of  successive  approximations 


u0(x)  =  Uf  (x), 
un+l(x)  =  H[un]  +  U  f(x) 


[Eq.  A-31] 
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will  converge  to  the  unique  solution  of  Equation  A-29.  This  sequence,  called  the 
Born  sequence,  has  been  used  to  discuss  effects  of  atmospheric  turbulence.  The 
approximation  of  the  solution  of  Equation  A-28  by  the  jth  term  in  the  sequence  of 
Equation  A-31  is  called  the  jth  Born  approximation. 

An  approximation  developed  by  Rytov  and  others  models  U(x)  as 

U(x)  =  Uf  (x)eVl+V2+-  [Eq.  A-32] 

where  vpj  is  taken  to  have  magnitude  0(81),  s  =  ||H||.  Substituting  the  right  side  of 
Equation  A-32  for  U  in  Equation  A-29  and  equating  terms  of  orders  s  and  s2  gives 

V,  =7TH[Uf]  [Eq.A-33] 

and 

[Eq.  A- 34] 

U  £  A 

Results  of  Prokhorov  [8]  and  others,  described  by  Hufnagel  [9],  show  that 
ensemble  averages  of  laser  beam  intensity  in  the  focal  plane  are  approximately 
Gaussian,  with  peak  (center)  intensity  <I>(0,  F0)  given  by2 

2 

<  I  >  (0,F0)  =  10- - - -y  [Eq.  A-35] 

fw[l  +  0.46Dw(2F0/kw)]3j 

where  10  is  the  intensity  at  the  center  of  the  initial  beam,  where 

w0=^,  w=  *  [Eq.A-36] 

'V  2  KW  q 

and  where 

1  _5  F  5 

Dw(2F0  /kw)  =  9.25k  jC2  (h(s))(F0 -s)3ds  .  [Eq.  A-37] 

s=0 

In  Equation  A-37,  s  is  arc  length  along  the  path  of  the  beam’s  axis,  h(s)  is  the 
height  of  the  axis,  and  C2  (h)  is  the  refractive  index  structure  function  of  the 

turbulent  atmosphere,  which,  as  the  notation  indicates,  we  assume  is  a  function  of 
height  only. 

2  Equation  A-35  corrects  a  typographical  error  in  Equation  35  of  LMI  report  [2], 
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The  beam’s  radial  intensity  profile  is  approximately  Gaussian,  with  e-folding 
length  Wturb  given  by 


wturb  =  w[l  +  0.46Dw(2F0  /kw)]s .  [Eq.  A-38] 

Using  methods  described  later  in  this  appendix,  we  evaluated  the  path  integrals  of 
Equation  A-37  numerically  to  calculate  the  effects  of  turbulence  on  a  laser  beam. 
We  also  used  the  method  in  the  section  on  extinction  to  generate  models  of  the 
atmospheric  extinction  coefficient  P  as  a  function  of  h  from  PcLnWin  outputs. 
Finally,  we  applied  the  blooming  calculations  of  Weichel  in  the  way  described  in 
the  section  on  blooming  to  determine  effects  of  blooming.  That  section  also 
includes  an  extended  discussion  of  the  turbulence-induced  variation  of  beam 
radius  with  distance  along  the  axis,  which  we  used  to  calculate  blooming. 

Extinction 


Our  basic  source  of  extinction  information  was  PcLnWin,  a  graphical  interface  to  the 
FASCODE  atmospheric  code  provided  by  Ontar  Corporation.  PcLnWin  is  fully 
documented  in  such  Ontar  publications  as  PcLnWin  Manual,  Version  3  (October 
1999). 

PcLnWin/FASCODE  requires  certain  user  inputs,  while  providing  default  values 
for  many  FASCODE  inputs.  We  used  the  following  inputs: 

♦  Line  file.  Using  data  provided  with  PcLnWin,  we  built  a  line  file  covering 
the  wavenumber  interval  from  7224.491  cm"1  to  7983.610  cm'1,  roughly 
±5  percent  from  the  chemical  oxygen  iodine  wavenumber  of  7604.523  cm' 

The  file  includes  N2,  O2,  H20,  C02,  O3,  N20,  CO,  and  CH4  molecules. 

♦  Atmosphere.  We  developed  results  for  three  of  PcLnWin’s  available 
atmosphere  models:  1976  Standard  Atmosphere,  and  two  models 
assembled  by  Ontar  from  various  sources  to  treat  the  winter  atmosphere  in 
mid-latitudes  (MLW),  and  the  summer  atmosphere  in  mid-latitudes 
(MLS). 

♦  Aerosol.  We  chose  FASCODE’s  Rural  aerosol.  This  choice  gives  23  km 
visibility  at  the  surface. 

♦  Troposphere/stratosphere.  This  transition  was  determined  by  the 
atmosphere  model  used. 

From  PcLnWin  outputs,  we  obtained  transmissivity  T  on  1-km  paths  at  heights 
from  0.1  to  15  km  above  the  surface.  These  data  allowed  us  to  tabulate  extinction 
coefficients  |3  at  the  observed  altitudes,  as 

P  =  -ln(T)km  1 . 
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We  then  obtained  extinction  coefficients  (3(h)  at  an  arbitrary  altitude  h  below 
1 5  km  by  linear  interpolation  in  the  logarithms  of  the  tabulated  extinction 
coefficients. 

Figure  A-4  compares  our  model  of  molecular  absorption  with  the  one  used  by 
Walter  and  Mani  [10]. 

Figure  A-4.  Comparing  Molecular  Absorption 


Figure  A-5  compares  our  extinction  coefficient  with  the  one  used  by  Walter  and 
Mani  [10]. 


Figure  A-5.  Extinction  Coefficient  Model  and  Data 


PcLnWin - Walter  a  Model 


A-12 


Calculations  of  Peak  Fluence  and  Target  Interaction 


Figures  A-4  and  A- 5  show  that,  while  PcLnWin’s  molecular  absorption  agrees 
closely  with  Walter’s  results,  PcLnWin’s  total  extinction,  and  thus  by  implication 
aerosol  absorption  and  extinction,  are  considerably  larger  for  heights  of  2  to  5  km 
than  those  of  the  model  used  by  Walter  and  Mani. 

Having  extinction  coefficients  as  functions  of  height  allowed  us  to  make 
numerical  calculations  of  the  transmissivity  along  a  path  of  length  L  on  which 
height  h(s)  is  given  as  a  function  of  path  length  s,  as  a  path  integral: 


We  evaluated  such  integrals  numerically,  using  the  method  described  below  in  the 
section  on  numerical  evaluation  of  path  integrals. 

Blooming 


The  beam  heats  the  air  through  which  it  passes,  decreasing  refractive  index 
unevenly  and  distorting  the  beam.  This  effect,  called  “blooming,”  reduces  fluence 
on  the  target.  Here  we  consider  this  effect. 


Citing  an  article  by  Weichel  [11],  the  APS  report  [6]  on  propagation  of  high- 
power  laser  beams  says  that,  for  a  Gaussian  beam  having  an  e-folding  radius  a(z) 
where  z  is  path  length,  the  ratio  of  the  intensity  at  the  target  with  blooming  to  that 
without  blooming  is 


I(Bloom)  _  1 

I(NoBloom)  1  +  (B0Ib)2 

where 

tj  _  dn0  a0PL 

Dq  —  ^ 

dT  47rp0w0ao 


and 


ib=7t1 


z=0 


a(z) 


1 


aow0q(z')T0 
o  a2  (z')w(z')a0T(z') 


-dz' 


dz 


[Eq.  A-39] 


[Eq.  A-40] 


.[Eq.  A-41] 


In  Equations  A-40  and  A-41,  a(z)  denotes  the  part  of  the  atmospheric  extinction 
coefficient  that  puts  power  into  the  air  through  which  the  beam  passes,  w(z)  the  effective 
wind  velocity  coming  from  the  motion  of  the  beam  relative  to  the  (possibly  moving) 
atmosphere,  T(z)  the  temperature,  and  a(z)  the  beam  radius,  all  varying  with  distance  z 
along  the  path.  The  subscript  0  on  any  of  these  denotes  the  corresponding  value  at  z  =  0. 
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We  used  the  model  of  atmospheric  extinction  described  above  in  the  section  on 

oc(z) 

extinction  for  the  normalized  extinction  coefficient  v  y .  The  following 


af 


a(z) 


subsections  explain  how  we  evaluated  normalized  beam  radius - -  and  wind 


nonnalized  wind  profile 


w(z) 


Wr 


Beam  Radius 


For  the  variation  of  beam  radius  a(z),  we  considered  the  variation  from  the  mirror 
radius  at  z  =  0  to  the  radius  of  the  turbulence-broadened,  focused  spot  at  z  =  L.  By 
“turbulence-broadened,  focused  spot”  we  mean  what  Fante  [12]  calls  the  “short¬ 
term  beam  spread”  for  the  focused  beam.  As  Fante  explains,  long-scale  turbulence 
causes  the  beam  to  wander,  while  short-scale  turbulence  spreads  it  out.  Thus  the 
target  receives  a  broadened,  wandering  spot  (Figure  A-6). 

Figure  A-6.  Broadened,  Wandering  Spot 


Source:  After  Fante  [12],  footnote  5. 


We  assume  that  the  heating  that  causes  blooming  is  influenced  by  the  broadening, 
but  not  by  the  wandering. 

Fante  [12]  gives  means  of  calculating  the  mean  e-folding  radius  of  the  broadened 
spot,  ps,  as  a  fraction  p(po)wull-b,  where 

3 

[Eq.  A-42] 


Po  = 


i 

1.46k2Ljc^L)(l-i;) 
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Calculations  of  Peak  Fluence  and  Target  Interaction 


When  po  is  either  much  less  than  or  much  greater  than  the  initial  beam  radius, 
Fante  [12]  provides  analytic  expressions  for  p(po).  As  shown  by  Figure  A-7, 
neither  of  those  conditions  is  likely  to  be  met  for  mirrors  of  the  size  likely  to  be 
used  in  laser  weapons  for  tactical  aircraft,  except  perhaps  at  very  low  altitudes  and 
extreme  ranges.  Rather,  in  the  situations  ALIRT  is  intended  to  treat,  po  is  roughly 
the  same  size  as  the  initial  beam  radius. 

Figure  A-7.  Variation  of  p0 


Fante  [12]  provides  numerical  means  of  calculating  ps  when  p0  is  roughly  the 
same  size  as  the  initial  beam  radius.  We  have  incorporated  Fante’s  method  for  this 
in  ALIRT,  by  making  a  numerical  model  of  Figure  3  of  [12]. 

Figure  A-8  compares  the  mean  radius  of  the  broadened,  wandering  beam  with  that 
of  the  broadened  beam  without  effects  of  wandering,  and  with  the  radius  of  the 
Gaussian  beam  propagating  from  the  laser  in  free  space,  for  a  horizontal  shot  from 
a  laser  weapon  at  7  km  height. 


3  Figure  A-7  also  shows  the  non-monotone  variation  of  Cn  with  height:  high-altitude  winds 
and  corresponding  turbulence  give  C“  a  maximum  at  around  10  km  height. 
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Figure  A-8.  Comparison  of  Beam  Radii,  Horizontal  Shot  at  7  km  Height 


— ♦ — Broadened,  wandering  Only  broadened  — * — Free-space 


The  figure  indicates  that  broadening  does  not  cause  substantial  departures  of 
short-time  average  beam  radius  from  the  free-space  case  at  heights  of  several 
kilometers  and  ranges  of  some  tens  of  kilometers,  though  wandering  does  cause 
the  long-time  average  beam  radius  to  be  significantly  larger  than  that  of  the  free- 
space  beam  in  such  cases. 

One  would  expect  greater  impact  of  turbulence  when  beams  propagate  through 
regions  of  greater  turbulence  intensity,  and  this  is  borne  out  by  calculations. 
Figure  A-9  shows  the  same  comparison  as  Figure  A-8,  but  for  a  horizontal  shot  at 
a  height  of  2  km.  The  figure  indicates  that  turbulence-induced  broadening  and 
wandering  are  significant  for  shots  at  lower  altitudes. 

Figure  A-9.  Comparison  of  Radii,  Horizontal  Shot  at  2  km  Height 
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Calculations  of  Peak  Fluence  and  Target  Interaction 


To  evaluate  the  iterated  integral  of  Equation  A-41  with  a  single  integral,  we  note 
that 


z=oa(z) 


aowoa(z')T0 

J=0a^(z')w(z')a0T(z') 


1 


-dz' 


dz=  J 


aoW0a(z')T0 


=o  a  (z')w(z')a0T(z') 


-dz 


dz' 


z=l  a(z) 

[Eq.  A-43] 


and  that  the  integral  with  respect  to  z  can  be  evaluated  numerically,  using  values 
of  the  variation  of  the  broadened  beam  along  the  beam  axis: 


j  — y-dz  =  g(z')  [Eq.A-44] 

z=L  a(z) 


With  the  g(z')  determined  in  Equation  A-44,  and  using  values  of  To/T(z')  from  the 
1976  Standard  Atmosphere,  we  evaluated  Ib  numerically  as 


Ib 


f  aoW0a(z')T„ 

J  2/  T/  l\  )  UZ- 

z'=o  a  (z  )w(z  )a0T(z  ) 


[Eq.  A-45] 


Wind  Profile 

We  needed  to  develop  a  wind  profile  w(z)  that  is  realistic,  and  for  which  the 
integral  of  Equation  A-45  converges,  and  for  which  Equation  D.29  of  APS  [6] 
will  apply.  That  equation,  specifically, 

5T 

pc  w —  =  al ,  [Eq.  A-46] 

t  dx 

is  the  starting  point  of  the  blooming  analysis  of  that  reference. 

For  air-to-air  engagements  (except  for  shots  made  flying  directly  toward  the  target 
or  away  from  it)  the  wind  relative  to  the  beam  will  be  large  along  the  entire  path, 
Equation  A-45  will  apply,  and  the  integral  will  certainly  converge. 

But  for  a  stationary  target  on  the  ground,  it  will  not  do  simply  to  take  the  velocity 
profile  generated  by  the  motion  of  the  beam  from  aircraft  to  a  fixed  ground  target. 
That  wind  profile  may  fall  linearly  to  zero  at  the  target,  and  the  integral  of 
Equation  A-38  would  not  converge. 

It  is,  however,  possible  that  typical  winds  would  be  large  enough  for  Equation  A- 
45  to  apply.  If  so,  then  a  simple  profile  w(z)  that  began  with  a  speed  like  that  of 
the  weapon-bearing  aircraft  at  z  =  0,  and  fell  to  a  typical  surface  wind  speed  at  z  = 
L,  could  be  used. 
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We  began  our  consideration  of  wind  profile  with  the  equation 


— +  q-VT  =  — V2T  +  — Q.  [Eq.A-47] 

dt  pcp  pcp 

expressing  conservation  of  energy  in  an  isobaric  motion  of  a  fluid  with  constant 
specific  heat  at  constant  pressure  cp  and  thennal  conductivity  k,  neglecting 
viscous  dissipation.  Equation  A-46  is  written  for  an  inertial  frame  moving  with 
the  beam  locally.  In  Equation  A-46,  q  is  fluid  velocity  in  cm/sec,  T  is  temperature 

in  °K,  Q  denotes  rate  of  heat  addition  in  watts/cm3,  k  is  thermal  conductivity  in 
watts/(cm  °K),  and  p  is  density  in  g/cnr . 

^  ^  q  /v 

Introducing  non-dimensional  variables  T,  q,  x  and  t  with  the  definitions 

-  _  kt 

Pcpd2 

q  =  q/w  [Eq.  A-48] 

x  =  x/d 
T  =  T  /  Tref 


where  x  is  position  coordinate,  w  is  fluid  velocity  relative  to  the  beam,  and  d  is 
beam  diameter,  brings  Equation  A-46  to 


5T 


wpcpd  3t 


+  q  •  VT  = 


wpc„d 


p  W 

V2T  +  - - Q. 


dT, 


ref 


[Eq.  A-49] 


In  Equation  A-49,  we  have  suppressed  the  A  symbols  for  clarity.  Equation  A-49 
shows  that  temperature  changes  will  satisfy 


q-VT  = 


w 

dTref 


Q 


[Eq.  A-50] 


provided  that - «  1 ,  or,  equivalently,  that  w  » - .  For  a 

wpCpd  pcpd 

representative  case,  k  =  2. 55x1  A-4  joule/(sec  cm  °K),  p  =  1.225xl0'3  g/cm3,  cp  = 
1  joule/(g  °K),  and  the  condition  that  Equation  A-49  governs  temperature 
changes,  and  w  »  0.01  cm/sec. 

This  condition  is  likely  to  be  met  by  ordinary  wind,  even  if  the  beam  is  fixed  in 
space,  as  it  would  be  on  a  stationary  target  on  the  ground. 
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Calculations  of  Peak  Fluence  and  Target  Interaction 


With  this  result  in  mind,  we  developed  the  profile  of  wind  relative  to  the  beam  in 
this  way:  If  VP  and  Vx  denote,  respectively,  the  velocity  vectors  of  platform  and 

target,  and  if  W  denotes  the  wind  velocity,  then  at  a  point  on  the  beam  located  a 
fraction  q  of  the  distance  from  platform  to  target,  the  beam  will  have  velocity 

V(g)  =  Vp+  g(VT  -VP)  +  W(g)) .  [Eq.  A-5 1] 

The  magnitude  Vperp  of  the  velocity  perpendicular  to  the  beam  will  be  given  by 
the  magnitude  of  the  cross  product  of  a  unit  vector  parallel  to  the  beam,  and  the 

beam’s  velocity  V  .  We  evaluated  that  cross  product,  arbitrarily  choosing  the 

components  of  W  to  be  a  constant  2  m/s  in  each  of  three  perpendicular  directions, 
to  obtain  w(z). 

In  addition  to  parameters  detennined  in  evaluating  IB,  evaluating  Bo  required  only 

a  value  for  —^2- .  We  used  the  value  -  lxlO'3  cm3/(g  °K)  given  in  the  APS 
Po  dT 
report  cited  above. 

Beam-Target  Interaction 


The  laser  beam  will  not  necessarily  strike  the  target’s  surface  perpendicularly,  and 
the  target  surface  will  not  necessarily  be  a  plane.  Here  we  consider  these  effects. 

For  perpendicular  incidence,  the  target  experiences  the  intensity  in  the  focal  plane 
X3  =  Fo.  Otherwise,  we  must  consider  variations  in  X3  as  well  as  in  xi  and  X2,  when 
X3  is  close  to  Fo. 

For  a  TEMoo  beam,  the  most  sensitive  variation  is  in  the  tenn 


1  1 

“1  ik“ 

— 3- -1 - A  — 3- - A 

W02  2  Wo  2 


where  A  =  1/Fo  -  I/X3  (this  follows  from  Equation  A-67  in  the  mathematical 
note). 


Now,  the  spot  size  will  limit  the  variation  of  X3  to  a  few  centimeters.  Accordingly, 

.  2kVWn2  „  2k2r2W2  .  „ 


the  term 


will  vary  from 


-  by  an  amount  of 


0(10 12  xlO  4xl0  2xl0  8  x  10  5),  which  will  cause  negligible  change  in  the 


quantity  of  Equation  A-52  .  Similar  considerations  show  that  the  tenn 
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k2Wp  A2  /  4  will  depart  negligibly  from  1,  while  X3  ranges  over  the  illuminated 
spot. 

Consequently,  over  the  part  of  the  target  surface  illuminated  by  the  beam,  the 
Poynting  vector  will  be 


P  = 


[Eq.  A-53] 


where  w  is  given  by  the  wturb  of  Equation  A-38  ande3  is  a  unit  vector  in  the  beam- 
propagation  direction. 

Let  us  first  consider  a  plane  target  surface.  Without  loss  of  generality,  for  a 
TEMoo  beam,  we  may  take  the  direction  of  xi  to  be  perpendicular  to  the  plane 
defined  by  the  beam-propagation  direction  and  the  normal  to  the  surface.  If,  then, 
yi  is  a  coordinate  in  the  plane  of  the  target,  parallel  to  xi,  and  yi  is  a  coordinate  in 
the  plane  of  the  target,  perpendicular  to  y  1,  in  the  target  plane,  the  intensity  will  be 


target  =COS0 - yexP 


7tW 


y2  +y2  cos2  9 
w2 


[Eq.  A-54] 


where  9  is  the  angle  between  the  normal  to  the  surface  and  the  beam-propagation 
direction.  Curves  of  constant  intensity  will  be  ellipses,  and  the  peak  intensity 
experienced  by  the  target  will  be  reduced  by  the  factor  cos9.  Figure  A- 19 
compares  1  kw/cm2  contours  at  normal  incidence  (circle)  and  at  69°  incidence 
(ellipse). 


Figure  A-10.  Comparison  of  Constant-Intensity  Contours 


If  the  target  surface  is  not  a  plane,  but  has  a  continuously  turning  tangent  over  the 
part  illuminated  by  the  beam,  then  it  may  be  parameterized  by 


y>  =Y,(yl,y2)  =  —°Sy  +  Y2Sin 


y)2  ,  (-yi 


R, 


■  + 


siny +  y2  cos 


y)2 


R. 


[Eq.  A-55] 
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Calculations  of  Peak  Fluence  and  Target  Interaction 


where  y3  is  a  coordinate  measured  perpendicular  to  the  plane  of  yi  and  y2,  where 
Ri  and  R2  are  the  principal  radii  of  curvature,  and  where  y  is  the  angle  between 
the  direction  of  yi  and  the  direction  of  the  plane  of  principal  curvature  Ri. 

With  Equation  A-54  and  Equation  A-55,  one  may  write  an  expression  for  the 
luminous  flux  into  the  surface  at  the  point  (yi,  y2,  Y-Tyi,  y2)).  The  result  will 
differ  from  Equation  A-52  for  constant  0  by  a  multiplicative  factor  whose 
difference  from  unity  is  bounded  in  magnitude  by  the  largest  value  of 

Vy?  +yl  m  the  illuminated  spot,  divided  by  the  smaller  of  Ri  and  R2,  and 
multiplied  by  the  largest  value  of  sin0. 

Since  Jyf  +  y\  is  typically  not  larger  than  about  10  cm  in  the  spot,  and  for 

reasonable  angles  of  incidence  sin0  is  0(  1 0"1 ),  the  multiplicative  factor  will  differ 
from  unity  by  O(10"2),  and  curvature  will  be  negligible,  if  the  smaller  principal 
radius  of  curvature  is  0(  1  m). 

It  seems  reasonable  to  assume  that  that  will  be  the  case,  when  the  target  is  a 
fighter  or  bomber  aircraft.  Accordingly,  we  will  not  consider  curvature  further 
here.  For  such  targets  as  air-to-air,  ground-to-air,  or  ground-to-ground  tactical 
missiles,  however,  curvature  effects  may  be  significant  and  should  be  considered. 

Numerical  Evaluation  of  Path  Integrals 


We  evaluated  path  integrals  such  as  the  one  in  Equation  A-37  numerically.  To  do 
so,  we  needed  to  compute  the  path  taken  by  the  laser’s  output  through  the 
atmosphere  undisturbed  by  turbulence. 

Effects  of  oxygen  and  nitrogen  molecules  dominate  the  atmosphere’s  refractive 
index.  Following  the  expression  for  n  given  by  Thomas  and  Duncan  [13],  we  used 

n  =  0.079P  /T  [Eq.  A-56] 

where  P  is  pressure  in  atmospheres  and  T  is  temperature  in  °K. 

For  MLS  and  MLW  atmospheres,  we  used  analytic  models  of  PcLnWin  refractive 
index  outputs.  Calculation  of  the  refracted  ray  requires  derivatives  of  refractive 
index  variation  with  height,  so  a  continuously  differentiable  model  is  desirable. 
We  made  the  models  by  linear  interpolation  in  the  divided  differences  of 
refractive  index  parameter,  and  integrated  the  resulting  piecewise  linear, 
continuous  function  with  an  added,  adjustable  constant  of  integration  to  obtain  a 
continuously  differentiable  model  of  refractive  index.  Figure  A-l  1  shows  the 
comparison  between  our  model  of  the  refractive  index  parameter  le6(n  -  1), 
where  n  is  refractive  index,  and  the  PcLnWin  data  for  MLS. 
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Figure  A-11.  Verification  of  MLS  Refractive  Index  Model 


Figure  A- 12  shows  the  comparison  between  divided  differences  of  the  refractive 
index  parameter,  and  the  derivative  of  our  model  of  the  refractive  index  parameter 
for  MLS.  The  tropopause  is  clearly  noticeable. 

Figure  A-12.  Verification  of  Derivative  of  MLS  Refractive  Index  Model 


— DD 
■  ■  Model 


The  refractive  index  variation  of  Thomas  and  Duncan  [13]  is  not  large;  it  is  about 
0.0003  near  the  surface.  Also,  tactical  laser  weapons  will  be  used  at  medium 
ranges,  generally  no  greater  than  100-200  km.  It  is  possible  that  using  straight 
paths  would  not  greatly  distort  integrals  like  that  of  Equation  A-37. 

Two  considerations  militate  against  neglecting  refraction,  however.  First, 
refraction  is  noticeable  at  extreme  ranges  and  low  weapon  altitudes.  Second, 

varies  sharply  with  height,  so  path  variations  may  have  significant  effects.  For 

example,  differs  by  more  than  10  percent  on  substantial  parts  of  the  two  paths 
shown  in  Figure  A- 13  (the  great  difference  between  the  scales  of  the  ordinate  and 
abscissa  of  the  figure  distorts  geometry). 
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Calculations  of  Peak  Fluence  and  Target  Interaction 


Figure  A-13.  Refracted  and  Straight  Paths 


—  Earth  - Straight  - Refracted 


Accordingly,  we  determined  refracted  paths  numerically.  We  did  so  by  integrating 
the  beam  equation  (from  Born  and  Wolf  [14]) 

—  (nr)  =  Vn 
ds 


and  the  relation 


dr 

- =  T 

ds 

between  the  vector  r(s) ,  which  gives  the  position  of  the  point  on  the  ray  at  path 
length  s  and  the  unit  tangent  vector  x  numerically,  using  the  difference  scheme 

ni+1fi+1  =  n1xi  +  (Vn)iAs 
Ii+1  =  I  ,As 


suggested  by  Berton  [15]. 

We  started  the  refracted  path  with  the  user-specified  altitude  and  elevation  angle, 
which  gave  usr(O)  andx(O) ,  respectively.  Please  note  that  starting  with  a  specified 
value  of  x  means  that  we  start  the  ray  in  the  direction  at  which  an  observer  looking 
with  the  laser  radiation  would  see  the  target,  not  the  direction  of  the  straight  line 
to  the  target  and  not  necessarily  the  direction  at  which  observer  looking  with 
radiation  of  a  different  wavelength  would  see  the  target. 


A-23 


Numerical  solution  of  the  beam  and  tangent  equations  gave  us  the  array  h(s,) 
specifying  beam  height  at  path  length  Si.  This  allowed  us  to  evaluate  the 
integrands  of  path  integrals  like  the  one  in  Equation  A-37,  which  have  the  form 
F(s,  h(s)),  and  then  to  approximate  the  integral  numerically  by  the  trapezoidal 
rule. 

Mathematical  Note 


To  evaluate 
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[Eq.  A-57] 


we  use  the  contour  integral  representation  for  Hermite  polynomials, 
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[Eq.  A-58] 


where  the  contour  C  is  a  circle  enclosing  z  =  0,  traversed  counter-clockwise 
(Abramowitz  and  Stegun  [16],  Equation  22.10.9).  The  representation  generates 
Hn(x)  for  arbitrary  complex  x. 

Substituting  the  contour  integral  in  Equation  A-58  for  Hn  in  Equation  A-57  and 
exchanging  the  order  of  integration  gives 
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[Eq.  A-59] 


where 
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[Eq.  A-60] 


The  inner  integral  of  Equation  A-59  is  readily  evaluated  by  completing  the 
square.  When  this  is  done,  we  find,  after  some  manipulation,  that 
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[Eq.  A-61] 


where 
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a  = 
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Considerable  further  manipulation  shows  that 
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and,  consequently,  that 
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Also, 
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so  that  the  argument  of  the  Hermite  polynomial  is  a  real  number. 
For  the  particular  case  n  =  0, 
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If  the  integral  is  evaluated  in  the  focal  plane  X3  =  Fo,  then  A  =  0,  and 
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[Eq.  A-62] 


[Eq.  A-63] 


[Eq.  A-64] 


[Eq.  A-65] 


[Eq.  A-66] 


[Eq.  A-67] 


[Eq.  A-68] 
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[Eq.  A-69] 


and 
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[Eq.  A-70] 


Using  the  result  of  Equation  A-70  to  evaluate  the  integral  expression  of  Equation 
A-21  for  U(xi,  x2,  F0)  when  Uo(xi,  x2)  is  given  by  Equation  A-22  leads,  after 
straightforward  manipulation,  to 
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[Eq.  A-71] 


where 


2F 

W  =  — °- .  [Eq.  A-72] 

kW0 

If,  finally,  the  sum  in  Equation  A-71  has  just  two  terms,  with  (poo  =  0,  cp22  =  jt,  and 
a  fraction  y  of  the  total  power  is  in  the  TEM22  mode,  then  the  focal  plane  intensity 
will  satisfy 


I  cc 


V2(i-Y)J-  J— n2(xi  /w)H2(x2  /w) 


exp 


[Eq.  A-73] 


where 
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[Eq.  A-74] 
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